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Recent measurements of the doping dependence of the London penetration depth X{x) at low T in 
' clean samples of isovalent BaFe2(Asi_a:Pi,)2 at T <C Tc [Hashimoto et al, Science 336, 1554 (2012)] 

' revealed a peak in A(x) near optimal doping x = 0.3. The observation of the peak at T ^ Tc, 

points to the existence of the quantum critical point (QCP) beneath the superconducting dome. We 
associate such a QCP with the onset of a spin-density-wave order and show that the renormalization 

Qof X{x) by critical magnetic fluctuations, gives rise to the observed feature. We argue that the case of 
pnictides is conceptually different from a one-component Gahlean invariant Fermi liquid, for which 
■ correlation effects do not cause the renormalization of the London penetration depth at T = 0. 

(N ■ 

I PACS numbers: 74.70.Xa, 74.40.Kb, 74.25.Bt, 74.25.Dw 

: 

O ' Introduction - Properties of iron-based superconductors (FeSCs) have been at the forefront of research activities 
^ in correlated electron community over the last few years [1-4]. These materials have multiple Fermi pockets with 
^ ' electron-like and hole-like dispersion of carriers. It is well established that superconductivity in FeSCs emerges in close 
proximity to a spin-density-wave (SDW) order, and superconducting (SC) critical temperature Tc has dome-shaped 
c/3 ' dependence on doping, with Tc maximum near the onset of SDW order [5-8]. 

Several groups [9] put forward the scenario that superconductivity in FeSCs has s"' symmetry and emerges 
because SDW fluctuations increase inter-pocket interaction, which is attractive for s"' gap symmetry, to a level 
when it overcomes intra-pocket repulsion. Likewise, SC fluctuations tend to increase the tendency towards SDW. 
Once the system develops long-range order, the situation changes because SDW and SC orders compete, and the 
, order which sets first tends to block the development of the other. For some system parameters, such a competition 
O ' gives rise to a homogeneous coexistence of SDW and SC orders in some range of dopings [10-12]. A homogeneous 
^ ] coexistence of SDW and SC orders has been detected in 122 materials - electron-doped Ba(Fei_xCo2;)2As2 [7, 8, 13- 
18] and hole-doped Bai_2:K2:Fe2As2 [19, 20]. In the third class of 122 materials - an isovalent BaFe2(Asi_2:Pa;)2, the 
coexistence between SDW and SC order has not been yet probed experimentally, but the odds are that the two orders 
K*" ' do coexists because the phase diagram of BaFe2(Asi_a;Pj:)2 is quite similar to that for Ba(Fei_a;Coa;)2As2 [23]. 

The coexistence implies that the SDW transition line extends into superconducting phase. If this line reaches T = 0, 
the system develops a magnetic quantum-critical point (QCP) beneath the superconducting dome [25], see Fig. 1. A 
' magnetic QCP without superconductivity has been analyzed in great detail [26, 27], and it is known that quantum 
fiuctuations near this point give rise to non-Fermi liquid (NFL) behavior and to singularities in various electronic 
characteristics. An SDW instability inside the d-wave SC state has been analyzed in [28] and was shown to give rise 
' to NFL behavior of nodal fermions. 
^-H ' Coexistence observed in FeSCs brings about the new issue of whether there are electronic singularities at a magnetic 
QCP which develops in the presence of an SC order. Of particular interest are the singularities in quantities such as 
the penetration depth A(a;), which measures electronic response averaged over the whole FS. Early experiments [29] 
on Ba(Fei_2;Coa;)2As2 found no special features in X{x) at the onset of SDW order, but recent measurements in 
^ BaFe2(Asi_^Pj;)2 (Ref. 24) found a peak in A(a;) at the smallest T <C Tc at around optimal doping (see inset of 
■ ■ Fig. 1). The authors of Ref. [24] speculated that the peak likely indicates that there is a QCP beneath a SC dome and 
argued that the peak in A(x) is a generic feature of 122 FeSCs, but it is more difficult to detect it in Ba(Fei_2:Coa;)As2 
because of greater degree of electronic disorder caused by Co doping. Another potential reason why the peak has 
been observed only in BaFe2(Asi_:rPx)2 is that this material possess gap nodes [30] what generally leads to stronger 
effects due to quantum fluctuations. 

In this communication, we analyze the behavior of X{x) under the assumption that the QCP is associated with the 
development of SDW order beneath a superconducting dome. A preemptive nematic order may also play a role [31], 
but we will not dwell into this. 

London penetration depth near QCP. - In general, the peak in A(x) at a SDW QCP can emerge by three reasons: 
(i) a non-monotonic behavior of A near a QCP already within the mean-field theory (like the peak in the specific heat 
jump at Tc at the onset of coexistence with SDW [32]); (ii) critical fluctuations at the onset of SDW, not specific to 
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the form of the gap; (iii) critical fluctuations specific to the presence of the gap nodes. Besides, A can either diverge 
at a QCP, or get enhanced but stay finite. It was found recently [33] that, at the mean field level, the variation of A 
is smooth and cannot explain sharp feature observed in Ref. [24]. Here we investigate the effects of critical magnetic 
fluctuations. We find that fiuctuations associated with SDW QCP beneath a SC dome give rise to the enhancement 
of the effective mass m* . The mass does not diverge because SC order cuts infrared singularities, but nevertheless 
m* /m at the QCP is noticeably enhanced. We argue that the enhancement of m* gives rise to a sharp peak in A(a:) 
at the onset of coexistence with SDW. We also find that the presence of the nodes in the gap is not sufficient to 
transform a peak into a divergence because the dominant contribution to m* comes from the region away from the 
nodes. 

London penetration depth in a type-II superconductor 
with cubic symmetry is expressed via zero-momentum com- 
ponent of the electromagnetic response tensor Qij{k) — 
{5ij — kikj/k'^)Q{k), which relates vector potential A and 
the current density j: ji{k) — —Qij{k)Aj{k). The temper- 
ature and doping dependent penetration depth is given by 
A~^(T, a:) = {Att / c)Q{T,x), where c is the velocity of light. 
The kernel Q{T,x) is related to current-current correlation 
function in the limit of zero frequency and vanishing momen- 
tum and is expressed via the superfiuid density (T, x) as 
Q = e^ns/mc, where m and e are the mass and the charge 
of an electron. Then A^ = rac^ / {Ane^ns). In Galilean in- 
variant, one-component fermionic system superfiuid density 
at T = is equal to the total density of fermions n{x). In 
this situation A(r = 0,x) does not depend on Fermi liquid 
corrections and remains the same as in a Fermi gas [34, 35]. 
Diagrammatically, superfiuid density is given by the sum of 
two bubble diagrams made out of normal and anomalous 
Green's function, and the independence of ^^(r — 0,x) on 
the electron-electron interaction is the result of the cancela- 
tion between self-energy and vertex corrections to these dia- 
grams. At T > the T— dependent part of Us does depend 
on Fermi liquid parameters [35]. 

We find, however, that in iron-pnictides the situation is 
different because these systems have multiple Fermi pockets, 
and s"' pairing originates from iter-pocket interaction. The 
interplay between self-energy and vertex corrections then de- 
pends on the orientation of Fermi velocities and the values 
of superconducting order parameters at different FSs. We 
find that self-energy and vertex corrections generally do not 
cancel, and the penetration depth is roughly proportional to 
m* /m. 

We followed earlier works [9] and assumed that the most 
relevant interaction in Fe-pnictides is between hole and electron pockets, separated by Q = (tt, tt) in the folded Brillouin 
zone, and that the gap has s"' symmetry and changes sign between electron and hole pockets. We calculated the 
leading interaction correction to A(a;) in the one-loop approximation. This perturbative analysis is justified because 
renormalized \{x) does not diverge even at an SDW QCP. There are 16 diagrams with one-loop corrections to 
current-current correlator, half of them are self-energy and half are vertex corrections. We evaluated the diagrams 
(see Supplementary material for details) and found that self-energy and vertex corrections are of the same order, and 
both decrease the superfiuid density and increase the penetration depth. To be brief, below we analyze how \{x) is 
affected by inserting fermionic self-energy into current correlation function. A straightforward calculation yields, at 
one-loop order 
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FIG. 1: [Color online] Lower panel: a theoretical phase 
diagram of 122-type iron-based superconductors in tem- 
perature va doping plane. Critical temperatures Ts and 
Tc indicate transitions into pure SDW and SC phases re- 
spectively. A QCP lies beneath of the superconducting 
dome and separates pure SC and coexistence SC-I-SDW 
phases. Upper panel: the theoretical behavior of the pen- 
etration depth A at T = 0. In the mean-field approxima- 
tion (dashed line), A diverges at the edges of the super- 
conducting dome, flat inside pure SC phase and increases 
monotonically as system moves towards the pure SDW 
phase. Beyond mean-field, magnetic fluctuations give rise 
to a peak in A at the onset of the SDW order (solid line). 
The peak in A has been observed in Ref. [24]. 



A2(r = 0,x) = A|cs 



(1) 
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where 



/?= ( V(l-Zj(fcF)) ) = (hma,„„S,(fcj^,u;„)\ . (2) 



Here j labels Fermi pockets, S is a diagonal (normal) self-energy, which generally depends on the location of kp on 
the corresponding FS, and (. . = J^^ . . . d(i)/2TT. In a situation when the dependence of T,j{k,ujrn) onk — kp can 
be neglected, the quasiparticlc residue is related to mass renormalization as Zj{kp) = 'm/m*(kp). 

A similar expression for A has been earlier obtained for heavy- fermion superconductor UBei3 [36], which is a two- 
component system of conduction d— electrons and localized /-electrons, of which only the first carry the current. It 
is tempting to extend the one-loop result (1) to oc ^^ m/m*, but we caution that non-cancellation of one-loop 
self-energy and vertex corrections to the current correlator docs not necessary imply that vertex corrections can be 
simply neglected. An example of more complex behavior beyond one- loop order has been recently considered in [37]. 

Evaluation of the ferrnionic self-energy. We consider the minimal three-band model of two elliptical electron FSs 
and one circular hole FS and with interactions (?sdw and gsc in the SC and SDW channels. The phase diagram of the 
model has been obtained before [12]. We focus on the region where at T = the system has a long-range SC order 
and is about to develop an SDW order. Renormalization of mass on all FSs is of the same order, and for brevity we 
show the calculations of m*/m for one FS. 

Potentially singular self-energy comes from the exchange of near-critical SDW fluctuations. In the normal state, 
these fluctuations are overdamped and are slow compared to electrons. In a SC state, the dynamical exponent changes 
from 2; = 2 to 2; = 1 because fermions which contribute to bosonic dynamics become massive particle-like excitations. 
Such systems have been earlier discussed in the context of cuprates [26] and we follow the same approach in deriving 
the expressions for the self-energy and spin polarization operator in the SC state in our case. 

The one-loop self-energy due to spin-fluctuation exchange is given by 

S,(fe,a;„) = / p^L{q,nrn)Qj{k -q,uJn- flm) (3) 

where Um = 27rT{n +1/2) and Qm = ^nrnT are ferrnionic and bosonic Matsubara frequencies respectively. The 
normal and anomalous components of the Green's function in the SC are 

a,(fe,a.„) = ^-;'2'2+l- 7-,(fe,c.„) = ^^^|^, (4) 

where = ^j(fc) = Vj^p{k — kp), and the energy gap A^- is equal to Ah, on the hole FS and Ae((/') = — Ae(l±acos20) 
on the two electrons FSs (we choose Ah, A,, > 0). The gaps on electron pockets have nodes when a> 1. 
The spin-fluctuation propagator is given by 

L{q, Clrn) = -1 , L ^ . . (5) 

where the polarization operator Il{q,Slrn) is 

n(,, n^) = N,TY:jd^^ ^ ^ ^ ^ . (6) 

Here oj± = a;„ ± flm/^, ^± = C i ^/2, and we replaced the integration over momentum fe by J . . .d^k/Aw"^ = 
Nf J . . . dS^d(j)/ {2tt), where Nf is the density of states. Parameter 5 = 5^ + 5q accounts for the doping-induced 
modification of the FSs. The term = 60 + 62 cos 20 describes changes in the FSs radii and overall shape (ellipticity) , 
while the term 5q = vpqcos{(l) — tp) describes the relative shift in the centers of FSs, where and ip are the directions 
of kp and q. The magnetic SDW critical point is determined in terms of doping parameters and 62 from the 
condition F = 0, where F = {g~^^ + U{0, 0))NJ^. 

We first consider the case of equal gaps on both FSs {a = 0, A/j = Ae = A) and then discuss how the results 
are modified in the case when the gaps on electron pockets have nodes. Earlier calculations show [11] that there is 
a broad parameter range 0.8 < S2/S0 < 4.7 for which SDW order emerges gradually, and its appearance does not 
destroy SC order, i.e. SDW and SC orders coexist over some range of dopings. Since we are interested in T = limit. 
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it is sufficient to evaluate the propagator of magnetic fluctuations Eq. (5) only at small frequencies and momenta. A 
straightforward expansion leads to (see Supplementary material) 



where 



r = In 



Tc,o 
Ts.o 



|<50|arccoshJl + ,52/A2 



(7) 



arccosh 



+ 



|<5^|(A2+<52)3/2 



(8b) 



and 



77(A,5,^) = -(cos2(</)-V') 



2A2_^2 arccosh ^/ 1 + ^2 /A 



(A2+52)2 



(8c) 



In Eq. (8a) we absorbed coupling constants ^sdw idsc) into the 
transitions into a pure SDW (SC) states. 

Without superconductivity, a magnetic transition at T = 
would be into an incommensurate phase, as evidenced by a 
negative sign of 77 for A = [12, 31]. In the presence of SC 
order, the commensurate (tt, tt) magnetic order is stabilized 
provided that relevant S'^ < A^, which we assume to hold. 
By order of magnitude, x ^ V ^ l/A^. 

Substituting Eqs. (4) and (7) into Eq. (2) and integrating 
explicitly over the momentum transfer q (see Supplementary 
material for details), we obtain the fermionic residue for a 
direction cj) along the Fermi surface in the form Z{(j)) = 1 — 
I{(j))F. Here /(0) accounts for the (non-singular) angular 
dependence and is normalized such that {(l)h) = 1, where (j)h 
is the direction of a hot spot {& kp point for which kp + Q 
is also on another FS, see Fig. 2(a), and F accounts for the 
dependence on the distance to hot spot, measured by F, and 
on the system parameters and 62- In explicit form F — 
(F(k(V')),7(V')))v„ where k(V) = x/vW, 7 = ^/vii^)^^ 
and 



(A2+52)5/2 

corresponding critical temperatures T^^ (Tc^o) for the 
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Using Eqs. (1) and (2) we find 



(9) 



(10) 
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FIG. 2: [Color online] Scaling function F(k, 7) which ac- 
counts for the interaction correction to the London pene- 
tration depth A'^/A^pg — 1 oc F(k,7) is plotted versus 7 
which measures the distance to the quantum critical point 
for three different combinations of the system parameters 
encoded by k = 0.1,0.25,0.5. (see text). Inserts - (a) 
hole (circular) and electron (elliptical) Fermi surfaces, (i>h 
marks the location of a hot spot; (b) the dependence of 
-F(k,0), normalized to the pre- factor in Eq. 9, on k. 



Because angular integrals over </) in /(0) and over '\\) va. F are non-singular, the dependence of /3 on the distance to 
the critical point and on system parameters can be approximated by /3 ^ i^(K, 7). It is apparent from the integral 
in (9) that F{k,j) is finite even in the limit 7—^0, which implies that the penetration depth remains finite at the 
SDW QCP. Still, F(k,7) is peaked at the SDW QCP (when 7 = 0), and decreases as F{k,j) cx ln7/^ at 7 > 1. 
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We illustrate this behavior in Fig. 2. Because S\ ^ F{K,'y), the penetration depth is also peaked at the QCP. This 
behavior is in agreement with the data for isovalent BaFe2(Asi_xPa;)2 [24]. 

By order of magnitude F{k,0) = 0(1), hence /3 = 0(1). The enhancement of = X%csi'^ + ^) at the SDW QCP 
is larger if magnetic order remains commensurate (tt, tt) even in the absence of superconductivity. In this situation, (5o 
and 82 are not restricted to be smaller than A, and, if they are larger, (3 is enhanced by ((5/A)^. We caution, however, 
that once j3 becomes large, the one-loop approximation is no longer applicable, and, in particular, vertex corrections 
has been analyzed in more detail [37]. 

We next computed F[n, 7) for the case when SC gap has nodes on electron FSs. We found that, roughly, the angular 
dependence of the gap renormalizes k downwards. This, however, does not change A qualitatively - at a magnetic 
QCP -F(k,0) increases when n decreases, but still remains finite. We illustrate this in Fig. 2(b). The reasoning is 
simple: the nodes of s+~ gap are located at accidental kp points which generally differ from hot spots. In the special 
case when the gap nodes coincide with hot spots, Z at a hot spot diverges logarithmically at a SDW QCP, but 
momentum integral of Z{(j)) is still finite, hence A remains finite even in this ease. 

Conclusions - In this paper we considered the behavior of the penetration depth A(a;) in a clean Fe-based s"' 
superconductor at the onset of a commensurate SDW order inside the SC phase at T = 0. We found that the 
penetration depth remains finite but has a peak at the onset of SDW order. The peak originates from the renormal- 
ization of the current-current correlator due to interaction with low-energy magnetic fluctuations. Both self-energy 
and vertex corrections give rise to a sharp increase of A near the QPC. The magnitude of the peak is larger when 
s"' gap has accidental nodes, but still remains finite at the onset of SDW order. Our results are in good agreement 
with the measurements of the penetration depth in the isovalent BaFe2(Asi_a;Pa;)2 inside the superconducting phase. 
Experiment [24] shows that that A has a peak at roughly the same doping where the Neel temperature Tg intersects 
with Tc- Although there are no direct measurements yet of SDW order inside the SC dome, our results strongly 
support the scenario that SDW order in BaFe2(Asi_a;Pa;)2 persists into SC phase, as it happens in other Fe-based 
superconductors, and that the peak in the penetration depth occurs at a magnetic quantum-critical point inside the 
SC dome. 

We thank Y. Matsuda, T. Shibauchi, R. Fernandes, S. Maiti, R. Prozorov, and S. Sachdev for useful discussions. 
A.L. acknowledge support from Michigan State University. M.G.V. is supported by NSF grant No. DMR 0955500. 
A.V.C. is supported by the DOE grant DE-FG02-ER46900. 



Derivation of F in Eq. (8a) 

Distance to the quantum critical point in our model is defined by F = (g"^^ - n(0,0))/7V/, where n(q, 17) is the 
polarization operator. Here we provide the details of the derivation of the equation for F [Eq. (8a) in the main text]. 
Prom Eq. (6) of the main text, we have: 

n(0,0) = -27rTiV, 5^ /-^\ , El = ^l + A\ (11) 



where (. . .)^ = j^^ . . . d(j)/2n. Hence 



The coupling constant gsdw can be eliminated in favor of the transition temperature into the pure SDW state, 
Tsfi = (2e'''^/7r)Aexp(— l/gsdwA^/)j which allows us to rewrite the equation above as 

Next we add and subtract 1/ En inside the angle brackets and take the zero temperature limit. The infrared convergent 
part of the Matsubara sum can be safely converted into the frequency integral, while infra-red divergent part is 
evaluated using 27rT [l/i^« - l/i^n] = ln(7re-'^^r/A). This yields 

r = m f / r 'I \ _ m f ^Ifl^) . (14) 
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Combining two logarithmic terms and using BCS expression between the gap and critical temperature A = TTe~'^'^Tcfi 
one finds for their difference hi(Tcfi/Tg^Q). The remaining energy integral can be computed analytically with the help 
of the formula 



/ 

Jo 



dx arccosh(6/a) 



(15) 



Combining last two expressions, we reproduce Eq. (8a) of the main text. 

Expansion of the polarization operator n(q, Q.rn) 

To obtain Eq. (7) of the main text, we need to expand the polarization operator n(qr, fi) to order and fi^. 
Expanding the integrand in Eq. (6) of the main text and integrating over ^ we obtain 

n(q, nm) = n(o, o) + (5n(o, Um) + 5ii{q, o), (i6) 

SU{q,0) = -N.TivMrj: J d<^^ iuSi + Sl + l^f ' ^''^ 

At T ^ 0, Matsubara sums over ujm can be converted into frequency integrals. Integrating over cOm, we reproduce 
Eqs. (7), (8b) and (8c) of the main text. 

Calculation of the self-energy diagram 

To calculate X/Xbcs, we need to know 1 — Zikp) = limt^^o diuj^T,j{kF,i^m)- The self-energy is introduced in Eq. 
(4) of the main text. We have 

i-z(..).3r.i.i:.(..o) ^g(%-;^;-") (18) 

Using = —din, taking the limit a; — )• 0, and integrating by parts, we obtain 

1 - Z{kF) = 3Tj2^^^S{k, - q, -n) (20) 

Let's first focus on fej? at a hot spot. Using the explicit forms of the SC Green's function from Eq. (4) and SDW 
fluctuation propagator from Eq. (7), we express the quasi-particle residue as 

1 - Zikp) = ^J2 (^^1^2 + + r)2(4g2 + Q2 + ^2) ^^^^ 

The position oi kp on the Fermi surface is specified by the angle (j), x and F are given by Eqs. (8a) and (8b) in the 
main text and depend on system parameters, and r] is given by Eqn. (8c) and depends on system parameters and on 
the direction of q set by angle (i.e., r] = r]{tjj)). Taking the limit of zero temperature and introducing dimensionless 
variables 

. = ^ = <*)-^y (-) 

we re-express Eq. (21) as 

l-ZikF) = {F{K{i;)Mm^ (23) 
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where 



i^(«,7) 



3k 



+00 



^dxdydz 



(a;2 +y2_^_ ^^2 + ^^(3.2 + ^2 + 1) ■ 



(24) 



Integrating over y and then over x, we obtain Eq. (9) of the main text. 

When kp is not at the hot spot, the expression for thc^ self-energy becomes more complex, because an intermediate 
fermion at fe^ + Q is no longer located at the Fermi suriuc x;. This reduces the self-energy compared to its value at a hot 
spot. Roughly, A in Eq. (21) get replaced by -y/A^ + vjpkj^{(f>), where kpicf)) is the deviation from a hot spot along the 
corresponding Fermi surface. With this modification, 1 — Z{kF) acquires an extra factor k. l/y'l + (kF{(t>)/koY, 
where fco = A/vf- At a hot spot, (f) = <ph and I{(l)h) = 1- The parameter 7 also has to be modified accordingly, 
but this is not important for our purposes because we treat 7 as a phenomenological parameter which measures the 
distance to a magnetic critical point. 

It is instructive to analyze in more detail F(k, 7) at the SDW QCP, when 7 = 0. We have 



F{k,0) = 



+CXD 



(1 - k)z^ + 1 



1 



arccosh ( ^f-i^} 



One can easily make sure that at k <C 1, the dominant contribution to the integral comes from small z ■ 
z, arccosh(-\/-z2 + l/y/K\z\) « \n{2/y/K\z\), and hence 



(25) 



K. For such 



F(k,0) 



With the logarithmic accuracy this yields 



/+00 
. -00 



Kz'^ In ^ 



^/k\z\ 



1 + ^2 



{Z^ + 1)2 



dz. 



K « 1. 



8TTr]Nfvj,A 

In the opposite limit k 1 the integral predominantly comes from z ~ 1/-\/K) and we have 

1 lUK 



F(k,0) 



UTT^TjNfvjrA v^' 



K > 1. 



(26) 



(27) 



(28) 



Interplay between self-energy and vertex corrections to the current correlation function 

In this section we present calciilations to prove the point in the main text that one loop self-energy and vertex 
corrections to current correlation function in a superconducting state of an Fe-pnictide add up rather than cancel. 
We argue that non-cancelation is the consequence of the fact that the system lacks Galilean invariance (electron 
pockets are located at the corners of the Brillouin zone), and it holds even if the dispersion within each band can 
be approximated by a quadratic one. We contrast Fe-pnictides to single-band superconductors with the quadratic 
dispersion. In the latter vertex and self-energy corrections cancel out, as required by Galilean invariance [34, 35, 38]. 
To emphasize this point, we assume that the electronic dispersion near hole and electron pockets is quadratic and 
consider two cases only intra-pocket interaction and only inter-pocket interaction. The first case models Galilean 
invariant case (hole and electron pockets do not couple), the second case models s+~ superconductivity in Fe-pnictides 
(hole and electron pockets are coupled, and the superconducting gaps on the two are determined self-consistently). 

We compute the interaction correction to the function Q{k) introduced in the main text. In the London limit, 
Q{k 0) = Q = -c^^ dr J dx{j^{T)jx{0)), where j = e X]p-,a V/'p,"' a = 1,2 accounts for spin 
projections. The penetration depth is related to Q as — {Ait/c)Q. 

In the non-interacting case, Q = Qqis given by the two bubble diagrams made of the normal and anomalous Green's 
function. Fig. 3. At T = we have 



8 



where f ± = C i the overall factor of 2 is due to spin summation, d is the inter-layer spacing, Np = m/(27r) is 

the 2D density of states, and the normal and anomalous Green's functions Q and IF are 



, IF{S,p, UJn) — 



A 



(30) 



.2+e|+A2' -^^P'-"^- ^2+^2+A2- 

We followed Ref. [39] and subtracted from Qo in (29) its expression in the normal state, (the second term, with 



(a) 



g 




(b) 



g 




FIG. 3: Diagrammatic representation of the current correlation function in a BCS superconductor. Diagrams (a) and (b) 
represent normal and anomalous contributions. 



fj„(^p-, cli„) — l/{iijjn — Cp)): because there is no superconducting current in the normal metal. In this form, the 
integrand is regular and the integrations over momentum and frequency can be performed in arbitrary order. We 
integrate (29) over ^ first. The second term then drops out, and in the first we can safely set fc = after the integration. 



The angular integration yields / (d6'/27r)w^ - — Vp/2 ~ p|,/2m^, and Qo becomes 



Qo 



duJr, 



A2 



{ujI + A2)3/2 



Integrating over Un we obtain 



Qo 



N, 



(31) 



(32) 



Using Np — m/(27r) and the relation between pp and the full density n — p'jp/{2Trd) we find (at T — 0) Ug ~ n and 
Qo = ne^ /{mc) - the textbook result. Alternatively, one could integrate over LOm first. Then the first term in Eq. (29) 
vanishes at T = 0, while the second one exactly reproduces Eq. (32) with n = Ug (see e.g., Ref. [40]). 

We now proceed with the interacting case. We label interaction correction to Qo as 5Q {Q = Qo + 5Q). The total 
SQ is the sum of self energy and vertex corrections, SQ = SQy. + 5Qv- At one- loop order, 5Q is represented by the 
three basic diagrams in Fig. 4. 




FIG. 4: Diagrammatic representation of the interaction corrections to the current correlation function. Diagrams (a) and (b) 
represent self energy corrections SQt., and (c) represents vertex correction SQv- The straight lines denote fermion propagators, 
which can be either normal or anomalous, and the wavy Une denotes the boson propagator either in the spin or in the density 
channel. 
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Galilean invariant case 



We start with the case when the fermion-fermion interaction is restricted to the intra-pocket one and consider a 
single pocket with a quadratic dispersion. We analyze separately the cases when the intra-pocket interaction is the 
density and the spin channels. 

To account for the overall factors of different diagrams with normal and anomalous Green's functions, it is convenient 

to treat interaction corrections to Q in the Nambu formalism. We introduce the Nambu spinor X^e„ ~ t^^L'-'^^^Ll ~ 

['^i't;p.«»>' We denote by Tx,Ty,Tz the Pauli matrices operating in Nambu space along with the unit matrix 

To = 1. The current operator in Nambu representation is j = eJ^pXf^p^oXp- The non-interacting Nambu Green 
function is 



G{P, OJn) 



(33) 



The Nambu Green function (33) is related to the normal and anomalous Green functions as 



G{p, LJn) 



(34) 



Interaction in the density channel 

Consider first the case when the dominant intra-pocket interaction is between fermion densities. Such an in- 
teraction is generally described by the effective action Sint = ~{^/'^)J2q,m^^n,„PQ,^mP-q,-^m^ where p{q,Clrn) = 
^I]p,c^„,a V4,i^„,aVv+q,c-'+f2™,a ( " = 1) 2 accounts for Spin). In the Nambu representation, = Y.p'x''p+g'TzXp- The 
function LP{q,Qrn) is an effective propagator of charge fluctuation which in the case of intra-pocket interaction is 
peaked at g = 0, = and decreases when each of these two parameters increases. 

The self energy correction to the static current correlation function is 



(35) 



where Tr stands for trace operation in Nambu space, and the factor of 2 accounts for two self-energy diagrams in 
Fig. 4. Once the trace over the Nambu indices is taken, dQ^ can be represented diagrammatically by graphs shown 
in the first two lines in Fig. 5. Because our primary goal is to understand when vertex and self-energy corrections 
cancel out and when they add up, we make an additional simplifying assumption that the characteristic momenta 
and frequency, above which charge propagator LP{q, Q^) decreases, are much smaller than A/vp and A, respectively. 
In this situation, the integrations over the fermion and boson energies and momenta separate, and Eq. (35) takes the 
form 

(5Qe = e^c-irf-^ig/s , (36) 

where 

LP=TY,L''iq,nm)=LP{x = 0,t = 0) (37) 

is an instantaneous interaction, and 

/s = 2T ^ Tr [u^,xG'(p, Wn)TzG{p, LOn)TzG{p, UJn)Vp,xG{p, w„) . (38) 

Taking the trace in Eq. (38) with Eq. (34) we obtain cxphcitly, 

/e = {4}T ^ [g^ + gt- 2T^] , (39) 

where Q-{p^ujn) = Gi—Pi —^n) a-nd frequency and momenta arguments are omitted for clarity. 
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FIG. 5: Diagrammatic representation of the self energy, Eq. (35), (first and second row) and vertex corrections, Eqs. (40), (41) 
(third and fourth rows). Fermion propagators with a single arrow stand for a normal Green functions, Q. The double arrowed 
lines designate the anomalous Green functions, J-. These contributions (with the proper signs) are obtained after taking the 
trace in the expressions for Aq in the Nambu formalism, Eqs. (35) and (41). 



For the vertex correction part we obtain, within the same approximation 



SQv = e^c-^d-^L^Iv , 



where 



UJn)vf,xG{p, Un)T^G{p, Un)Vj;^xG{p, Un) 



(40) 



(41) 



Taking again the trace in Nambu space, we obtain the diagrams presented in the last two lines in Fig. 5. Evaluating 
them, we find 



(42) 



At T = the frequency summation is replaced by integration. Evaluating frequency and momentum integrals by, 
e.g., by using polar coordinates, we obtain 

Nf 



6A2 



r^g4^r^gl = r^ j-^g^ = tJ2 = « • (43) 



This gives 



ly. = -Iv = - 



1 



3 A2 



(44) 
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We see that sclf-cncrgy and vertex corrections to the static current correlation function in an s-wave superconductor 
cancel each other, as it indeed should be the case because of Galilean invariance. 



Interaction in the spin channel 

We now check that the cancelation in the Galilean invariant case holds also for the case when the intra-pocket 
interaction in the spin channel. The effective action is now Sint = ~{^/'^)J2q,m^'^,n,n'^Q,^m'^-q,-^m- The Nambu 
formalism is inconvenient in this case as the spin operator doesn't have a simple form of a product of a Nambu field 
operator and its conjugate. This formal difficulty can be alleviated by introducing Balian-Werthammer (BW) [41] 
extension of the Nambu formalism. The four-dimensional BW spinors are defined as 

= [*ii, *i2, *2i, *22]*'' = [cp-t, Cp-4, -clp^, . (45) 

These spinors form a tensor product of Nambu space and original spin space. Let ax,a-y,az be spin Pauli matrices, 
and (To the unit matrix. The current, density, and spin density operators are p = (1/2) X^^^'tp^'^, j = ev, v = 

(1/2) and a- = (1/2) Ep where 

Pp = Tz Cgi (To , Vp= vfTo 'Sxtq, a = To'Sicr , (46) 

and (g) stands for a tensor product. The Green function in BW space is expressed via the Nambu Green function, Eq. 
(33), as 

G = G (g) (TO . (47) 

The main advantage of the BW formalism is that it allows one to use standard diagrammatic rules, including (— 1) 
rule for each closed loop. 

We found that Eqs (36), (37), and (40) preserve their form, with L"^ instead of L^, but Eqs. (38) and (41) become, 
respectively 

/e = T ^ Tr [Vp,xG{p, 0Jn)&G{p, uJn)^G{p, uJn)%xG{p, ujn)] (48) 

and 

Iv = \tY, Tr [G(p, uJn)Vp,xG{p, LO„)&G{p, w„)F^,^G(p, a;„)<T] . (49) 

An extra factor of 1/2 as compared to (38) and (41) compensates for the double counting of degrees of freedom in 
BW formalism. Taking the trace in (48) and (49) we obtain 

7s = 3(tT|)T J2 [2-^* + {^^ -GG-+gl)+G^ + Gt] , (50) 



ly = Z^Lt [2^4 ^ 4j^2 I^g2 _ gg_ ^ g2_^ ^ gi ^ ^4 j ^^-^^ 

respectively. Using Eq. (43) we immediately find that 

I^=Iv = 0, (52) 
i.e., self energy and vertex corrections simply vanish. 
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The Ccise of inter-pocket interaction 



The basic formalism set up in the previous subsection can be applied to the case when the interactions LP''^{q, flm) 
are peaked at large momentum q = Q, and couple hole and electron pockets. Wc consider the simplest situation 
of perfect nesting, when the electron and the hole dispersions are = — ^^^^ = £,p. In this case, the velocity of a 
fermions at a point k on the hole pocket is opposite to that of a fermion at k + Q on the electron pocket. We also 
assume that Ag = —Ah = A. Finally, like in the previous case, we assume that LP''^{q,flrn) are rapidly decaying 
functions oi q — Q and flm, and that characteristic q — Q and flm much smaller than A/vf and A, respectively. 
This again allows us to factorize the integrations over fermion energies and momenta and over boson energy and 
momentum. 

Interaction in the density channel 
The electron and hole Nambu Green functions are 

Equation (36) should be modified to include both hole and electron pockets and becomes 

5gi: = ^i^(7£+7^), (54) 



where 



and 



LP=T J2 L^'iqMm) (55) 

qRiQ,Qm 



je^n, = 2r ^ Tr [vf;:}G<^\p,oj^)T,G'^^^\p,oj,,)T,G<'^\v^^^^ . (56) 



Taking the trace and counting pre-factors, we obtain 



= = (tr|)r ^ [-2j-3j-ft + j^^{2gegh - QeQh,- - Qe,-Qh + 2ge,-gh,-) 

-2:F^Fh [Gl - GeGe,- + Gl_) + GlGh + G%_Gh,^ ■ (57) 



The non-zero contributions are 

T Yl ^eJ'h = '2TJ2 ^iQeGh = 2T ^ TlGe,-Gh,- = 2T ^ TeJ'hGeGe,- = " ^ • (58) 



P,i^n P,W„ P,Wn P,W„ 

Substituting this into (57), we obtain 



J| = = Py^ = 4^ = -g^^- (59) 



The vertex correction contribution so the current-current correlation function is 



where now 

jeh{he) 



5Qv = ^Kilv + Iv) , (60) 



jen,ne> = T ^ Tr [G<^\p,UJr.)v^'J G<'^\p,OJ^)t,G'^^^\p,^^^^ . (61) 
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After taking the trace, this becomes 

2 



(vi) 

jeh _ jhe _ _} P' , 



T ^ [glGl + gl_gl_ + 2Tin + ^{Ql + Ql-) + + Ql-) 

In 

- IThT^ {GeGh + Ge-Qh- - Qe-Gh - QeGh-) ■ (62) 

The non-zero terms in (62) are 

T J2 SeGh =TJ2 Se,-Gl- =TJ2 ^eH = 2T ^ J^hJ^eGeGh = 2T ^ J^hJ'eGe,-Gh,- = ^ • (63) 



Substituting into (62), we obtain 



Evidently, SQ^ = 5Qv, i-c, sclf-cncrgy and vertex corrections add up rather than cancel out. By obvious reasons, 
non-cancelation survives even if we relax the assumption of a perfect nesting and of exactly opposite values of the 
gaps on hole and electron pockets, and do not require that L''''^{q,flrn) are rapidly decaying functions oi q — Q and 

■ 

The total 5Q = 5Qy, + SQv is 

The relative correction to the penetration depth 

6X SQ LP 



Ao 2Qo 6A2 ' 



(66) 



where Qo is given by Eq. (32), and the factor of 2 is because we include both electron and hole pockets. We see 
that the corrections tend to increase the penetration depth. The instantaneous interaction LP, given by Eq. (37), 
is of order Hinitrei, where Hint is the typical local interaction strength, and tret is the retardation time. Equation 
(66) shows that the perturbation theory is valid provided. Hint ^{tret^)- Because we assumed tret ^ the 
perturbative approach is valid even if the interaction exceeds A. 

Interaction in the spin channel 

Consider next inter-pocket interaction in the spin channel. Applying the BW formalism, we obtain, after straight- 
forward calculations 

je^W = T ^ Tr [v;l^'^G<'^\p,u;n)&G'^^'\p,u;n)&G<'^^ , (67) 

and 

je^hihe) = It ^ Tr [G<'^\p,UJr.K^^^^^G<'^\p,OJn)&&'^^^^ ■ (68) 

Taking the trace over spin and Nambu indices we get 

11: = 4= 3(4)T ^ [2J-3j-^ + 2J^eMGe - GeGe,- + GI,-) + GIGh + GI-Gh,- 

+ J'li^GeGh - Ge,-Gh - GeGh,- + 2Ge,-Gh,-)\ , (69) 
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I^h = Jhe ^ _3 ^ \Tl{gl + gl_) + glGl + 2TeTh{QeQh + Qe,-Qh,- - Qe,-Qh - QeGh,-) 



(70) 



The non-zero terms in Eqs. (69) and (70) axe listed in Eqs. (58) and (63). Substituting them, we obtain 



je jh jeh jhe 



(71) 



2A2 



We see that self-energy and vertex corrections again add up. Combining the two, we obtain 



5Q = 




(72) 



The correction to the penetration depth is 



5\ 
Ao 



5Q 
2Qo 



2A2 ■ 



(73) 



Aslamazov-Larkin diagrams 



In the calculations of the current correlation function in the normal state, it is important to include into consideration 
also Aslamazov-Larkin diagrams [42], which contain two boson propagators L coupled via two triangular fermionic 
vertex blocks. We analyzed these diagrams for our case and found that they are irrelevant. First, for the case when 
the gaps on hole and electron FSs are of equal magnitude (and opposite sign), there is a cancelation of the dominant 
Aslamazov-Larkin terms. The cancelation occurs at the level of the valuation of the fermionic triangular vertex - for 
each block there are two ways to arrange electron and hole Green's function lines and their corresponding momenta 
which cancel each other. For non-equal gap magnitudes, the cancelation is not prefect. Still, even in this case, the 
Aslamazov-Larkin terms can be safely neglected. The reason is the following: the overall factor in Aslamazov-Larkin 
terms contains two additional powers of the coupling constant. In the normal state, these two additional powers and 
eliminated by the need to regularize the divergence coming out of integration of the two bosonic propagators. In the 
superconducting state, the corresponding integral is finite because for small q each triangular block scales as Vpq/^- 
As a result, the two extra powers of the coupling do not cancel. 
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